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Inclusive production via gluon fusion will be the most important production channel for Higgs 
boson discovery at the LHC. I report on a calculation of the inclusive Higgs boson production 
cross section at pp and pp colliders at next-to-next-to- leading order in QCD. 



1 Introduction 

The Standard Model is almost thirty-five years old, and its basic assertion, that the weak and 
electromagnetic interactions unify in a spontaneously broken SU(2)l (g) U(1)y gauge theory has 
been spectacularly confirmed at the quantum level by the precision experiments at LEP and 
the SLC. Still, the crucial ingredient, the agent of electroweak symmetry breaking, remains a 
mystery. 

The benchmark for studies of the symmetry breaking sector is the minimal standard model, 
in which a single complex scalar doublet is introduced. Upon, spontaneous symmetry breakdown, 
a single Higgs boson is left as the signature of the symmetry breaking sector. The final run at 
LEP established a 95% confidence-level lower limit for the Higgs boson of 114.1 GeVEJ. The Higgs 
boson can also be constrained by its effect on precisely measured electroweak observables. The 
current best fit produced by the LEP Electroweak Working Group predicts that Mh = 85^^ 
GeV, with a 95% confidence-level upper limit of Mh < 196 GeV. 

With the end of LEP, the search for the Higgs boson is left to hadron colliders, the Fermilab 
Tevatron and the CERN LHC. Higgs production at hadron colliders is dominated by the gluon 
fusion mechanism, where gluons excite virtual top quark loops which couple strongly to the 
Higgs. Despite its dominance, this production process cannot be exploited at the Tevatron 
unless the Higgs boson mass is near the WW threshold. For lighter Higgs masses, the dominant 
decay mode, H — ► bb is overwhelmed by the QCD background and the overall rate is too low to 



permit the use of rare decay modes. Therefore, the Tevatron experiments must primarily rely 
on the associated production mechanism, qq — > H + W, for their Higgs search. 

At the LHC, however, gluon fusion is the most important production mechanism for Higgs 
discovery for masses below ~ 700 GeV, since the rate will be sufficiently high that one can use 
the rare decay mode H — > 77 in the low mass region. The discovery of the Higgs boson and 
the subsequent study of its properties therefore relies on a solid theoretical understanding of the 
gluon fusion production mechanism. 

Unfortunately, next-to-leading order (NLO) studies Btl of inclusive Higgs production do not 
provide this solid understanding. The NLO corrections are very large, of order 70 — 100%, 
and initial estimates of next-to-next-to-leading order (NNLO) corrections □ based on soft plus 
collinear resummation indicated that they too might be very large. The unsettled nature of such 
an important signal clearly calls for a renewed effort to bring this process under control. 

Last year, two groups presented calculations of the soft plus virtual contributions to the 
NNLO correction iH. It was known from NLO corrections that soft plus virtual terms alone 
significantly underestimate the full correction. The leading hard scattering term is completely 
collinear in nature and can be reliably obtained from resummation. The combined soft plus 
virtual plus collinear (SVC) approximation £j predicts that the total NNLO correction will be 
much more moderate than the initial estimate @. Still, one would like a full NNLO calculation 
to verify that inclusive Higgs boson production is under perturbative control. In this talk, I will 
present the results of that calculation 13. 



2 The Calculation 



In the limit that all quark masses except that of the top quark vanish, gluons couple to Higgs 
only via top quark loops. This coupling can be approximated by an effective Lagrangian Q 
corresponding to the limit mt — > 00, which is valid for a large range of Mjj, including the 
currently favored region between 100 and 200 GeV. The effective Lagrangian is 

£eff = "f C ^)G%G a ^, (1) 

where G" is the gluon field strength tensor, H is the Higgs field, v ~ 246 GeV is the vacuum 
expectation value of theHiggs field and C\(a s ) is the Wilson coefficient, which for this calculation 
we need to order (af)Ha. 

Three classes of Feynman diagrams must be evaluated to compute the NNLO cross section: 
(i) two-loop virtual diagrams for gg — > H; (ii) one-loop single real emission diagrams for gg — > 
Hg, gq — > Hq, and qq — > Hg; (hi) tree-level double real emission diagrams for gg — > Hgg, 
gg — > Hqq, gq — > Hgq, qq — » Hqq, qq — > Hgg, and qq — » Hqq. Sample diagrams are shown in 
Figure [j]. 
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Figure 1: Sample diagrams of (a) virtual (b) single real emission and (c,d) double real emission corrections. 



The most difficult part of the calculation, by far, is the computation of the corrections 
due to double real emission. While the (tree-level) matrix elements are quite easy to compute, 
the integration over double-emission phase space is quite complicated. We find, however, that 



if we expand the double real emission integrals as power series about soft limit (where the 



partonic center-of-mass energy is close to the Higgs mass, Mfj/s 



1), the integrals simplify 



dramatically. Moreover, the series expansion is well-behaved and converges quite rapidly 

Thus, while we have computed the virtual corrections, single real emission and the effects of 
mass factorization in closed analytic form, we have computed double real emission in the form 
of a power series. The closed-form results can also be readily expressed as power series, so we 
present the full partonic cross section as an expansion in (1 — x) and ln(l — x): 
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Note that if all coefficients are computed, this is an exact expression for the partonic cross section. 
In practice, we compute only a finite number of terms. The soft plus virtual approximation 
includes only the and terms at second order. The SVC approximation also includes the 
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coefficient. We have now computed all coefficients cf? through I = 160. As can be seen in 



Figure |3|, this is more than enough terms to obtain a reliable result for the total cross section. 

In Figure §(a), we show the cross section at LO, NLO and NNLO. At each order, we use 
the corresponding MRST parton distribution set0. One immediately sees that the true NNLO 
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Figure 2: (a) LO (dotted), NLO (dashed) and NNLO (solid) cross sections for Higgs production at the 
LHC (fiF = 1-tR = Mb)- In each case, we weight the cross section by the ratio of the LO cross 
section in the full theory (Mt = 175 GeV) to the LO cross section in the effective theory (Eq. (^)). 
(b) Scale dependence at the LHC. The lower curve of each pair corresponds to fiR = 2Mh, = Mij/2, the upper 
to hr = Mh /2, /if = 2Mh- The iC-factor is computed with respect to the LO cross section at fiR = /if = Mh- 



correction, while substantial, is much smaller than the NLO correction. Indeed, it is even a bit 
smaller than predicted by the SVC approximation. Nonetheless, it verifies that the SVC is a 
good approximation of the total cross section. 

Figure |2|(b) shows the renormalization and factorization scale dependence of the "K factor" , 
the ratio of the NLO and NNLO cross sections to the leading order cross section. The scale 
dependence of the NNLO cross section is still quite large, though somewhat smaller than at 
NLO. 

Figure H shows the rapid convergence of the power series expansion in (1 — x). Observe that 
the purely soft contributions underestimate the cross section by ~ 10% — 15%, while the next 
term, oc (1 — x) , overestimates it by about 5%. By the time the third term in the series is 
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Figure 3: 7f-factor for Higgs production at the LHC. Each line corresponds to a different order in the expansion 
in (1 — x). The renormalization and factorization scales are set to Mh- 



included (oc (1 — x) 1 ), one is within 1% of the result obtained by computing the first 18 terms 
(through (1 — x) 16 ). In light of the large scale uncertainty in the result, there is little or no 
precision to be gained from computing higher order terms in the expansion. 

In conclusion, we have computed the full NNLO corrections to inclusive Higgs boson pro- 
duction at hadron colliders. We find reasonable perturbative convergence and reduced scale 
dependence. 
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